Phase diagram of spin 1 antiferromagnetic Bose-Einstein condensates 
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We study experimentally the equilibrium phase diagram of a spin 1 Bose-Einstein condensate with 
antiferromagnetic interactions, in a regime where spin and spatial degrees of freedom are decoupled. 
For a given total magnetization m z , we observe for low magnetic fields an "antiferromagnetic" phase 
where atoms condense in the m — ±1 Zeeman states, and occupation of the m = state is sup- 
pressed. Conversely, for large enough magnetic fields, a phase transition to a "broken axisymmetry" 
phase takes place: The m = component becomes populated and rises sharply above a critical field 
B c (m z ). This behavior results from the competition between antiferromagnetic spin-dependent in- 
teractions (dominant at low fields) and the quadratic Zeeman energy (dominant at large fields). We 
compare the measured B c as well as the global shape of the phase diagram with mean-field theory, 
and find good quantitative agreement. 

PACS numbers: 03.75.-b,03.75.Hh 



One of the most active topics in the field of ultra cold 
quantum gases is the study of interacting many-body 
systems with spin Atoms with arbitrary Zeeman 

structure can be trapped using far-detuned optical traps. 
Quantum gases of fermions with spin larger than 1/2 



17] 



[4 1 and bosons with spin 1 (grill, 2 [g, llq l. or 3 
have been demonstrated experimentally. This opens a 
whole class of new experiments with spinful many-body 
systems, such as observation of squeezing among the dif- 
ferent spin components, coherent spin mixin g dy namics 
analogous to an internal Josephson effect [7-12, 16 1, or 
the study of sudden quenches across magnetic phase tran- 
sitions [3, E3]. 

The simplest example is the spin-1 Bose gas, where the 
spin-dependent interactions can favor either ferromag- 
netic (the case of atomic 87 Rb 0) or antiferromagnetic 
(the case of atomic 23 Na 0]) behavior, leading to different 
equilibrium phases. An additional but essential feature in 
experiments with gases of alkali atoms is the conservation 
of the longitudinal magnetization m z = n+i — n_i; Here 
n m denotes the relative populations of the Zeeman state 
labeled by the magnetic quantum number m = 0, ±1. 
This follows from the spin rotational symmetry of short- 
ranged two-body interactions : The only possible spin- 
changing two-body process is 



m = + m = — > m = +1 + m = — 1, 



(1) 



where two m = atoms collide to yield one atom in each 
state m = ±1 (or vice- versa), leaving m z unchanged. 
For antiferromagnetic interactions, the spin energy of the 
right hand side in Eq. (TT]) is lower than that of the left 
hand side. In most physical systems, the magnetization 
would relax by coupling to an external environment. In 
contrast, quantum gases are almost perfectly isolated and 
the conservation of magnetization plays a major role [3l| . 

In spite of intense theoretical activity Q, the equilib- 
rium properties of spinor gases remain relatively unex- 



plored experimentally. Most experimental work so far 
have focused on dynamical properties. For ferromag- 
netic Rubidium condensates, a recent experimental study 
concluded that the time needed to reach an equilibrium 
state, typically several seconds or tens of seconds, could 
easily exceed the condensate lifetime [l5|. For antiferro- 
magnetic 23 Na, the stationary regime after damping of 
spin-mixing oscillations has been studied for relatively 
high magnetization (m z < 0.5) 12|. Here also, long equi- 
libration times on the order of 10 s were observed. Both 
experiments worked with condensates with large atom 
numbers, well in the Thomas-Fermi regime, where spin 
domains are expected and observed in transient regimes. 

In this Letter, we present an experimental study of the 
phase diagram of spin 1 Sodium Bose-Einstein conden- 
sates with antiferromagnetic interactions. We work with 
small atomic samples containing a few thousands atoms 
held in a tightly focused optical trap; In this regime, spin 
domains are energetically costly, and spatial and spin de- 
grees of freedom are largely decoupled. We prepare the 
sample well above the condensation temperature with a 
well-defined longitudinal magnetization and no spin co- 
herence. At the end of the cooling stage, equilibration 
times of 3 s are used to ensure that thermal equilibrium 
is reached. We find, in agreement with theoretical pre- 
dictions, a phase transition from an "antiferromagnetic" 
phase where only the m = ±1 Zeeman components are 
populated to a mixed "broken axisymmetry" phase where 
all three Zeeman states can coexist. We determine the 
phase boundary and the shape of the phase diagram ver- 
sus applied magnetic field and magnetization by measur- 
ing the population of the m — state. Our measurements 
can be explained quantitatively by mean-field theory in 
the single-mode regime, where the atoms condense in the 
same spatial wave function irrespective of their internal 
state. 

We work with Sodium atoms cooled deeply in the 
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FIG. 1: (Color online) a: Absorption image of a spin 1 
BEC after expansion in a magnetic gradient, b: Horizontal 
cuts through the images in (a). The same function (shown by 
straight lines), only recentered and reweighted, is used to fit 
the density profile of each Zeeman state. 



quantum deg enerate regime using an all-optical cooling 
sequence li[ 1|| . In order to prepare the sample with a 
well-defined longitudinal magnetization and no spin co- 
herences, we start from a cold cloud in a crossed optical 
dipole trap loaded from a magneto-optical trap , with 
a magnetization m z w 0.6 resulting from the laser cooling 
process. To obtain higher degrees of spin polarization, we 
perform evaporative cooling in the presence of a vertical 
magnetic field gradient for about 1 s. Each Zeeman state 
sees a slightly different potential depth. Because of the 
combined action of gravity and of the magnetic gradient, 
evaporative cooling in this configuration favors the Zee- 
man state with the higher trap depth [20]. This results 
in partially or almost fully polarized samples with mag- 
netization up to m z 0.85. To obtain lesser degrees of 
polarization than the initial value m z « 0.6, we remove 
the gradient and apply instead an additional oscillating 
field resonant at the Larmor frequency. The two proce- 
dures together allow to prepare well-defined magnetiza- 
tions ranging from to ~ 0.85 with good reproducibil- 
ity and keeping the same evaporative cooling ramp in all 
cases. After this "spin distillation" , we transfer the cloud 
in the final crossed dipole trap and resume evaporative 
cooling (see the Supplemental Informations for more de- 
tails about the procedure). 

After the evaporation ramp, we obtain quasi-pure spin 
1 Bose-Einstein condensates (BEC) containing N ps 5000 
atoms in a trap with average frequency uj « 2ir x 0.7 kHz. 
To ensure that the cloud has reached a steady state, we 
allow for an additional hold time of 3 s after the evapo- 
ration ramp. We have investigated the dynamics of the 
spin populations as this hold time is varied for several 
values of magnetization and applied magnetic field. We 
found that the populations relaxed to steady-state values 
with a characteristic (1/e) time smaller than 1 s, much 
less than the finite lifetime of our sample, around 10 s. 

The populations of the Zeeman states m = 0, ±1 are 



analyzed after expansion in a magnetic field gradient pro- 
ducing a Stern- Gerlach force that accelerates atoms in 
m = ±1 in opposite directions. After a given expansion 
time (typically t « 3.5 ms), we take an absorption picture 
of the clouds (see Fig.QJi), and count the normalized pop- 
ulations n m of the Zeeman state to = 0, ±1. Note that 
the condensate is in a regime intermediate between the 
ideal gas and the Thomas-Fermi limits (with a chemical 
potential /i ~ 4fuU). 

For a Bose-Einstein condensate held in a tight trap 
as in our experiment, the energetic cost of spin domains 
is large (comparable to fuv per atom, much larger than 
the spin-dependent interaction energy). In this limit, 
it is reasonable to make the single mode approximation 

which 



(SMA) for the condensate wavefunction [2 lL |22j. 
amounts to consider that all atoms share the same spatial 
wavefunction independently of their internal state; The 
condensate spin remains as degree of freedom. To sup- 
port this approximation, we note that absorption images 
as in Fig. [TJt do not reveal any spatial structures or spin 
domains. Furthermore, we compare in Fig. [TJd the ob- 
served distributions with a common mode distribution. 
This common mode function is extracted from a Gaus- 
sian fit to the more populated cloud (m = +1 in this 
example), and then recentered and reweighed to match 
the populations of the other Zeeman states. We find very 
good agreement between the three spatial distributions 
in the whole range of parameters explored, and conclude 
that the SMA is indeed a good approximation in our case. 

Because the longitudinal magnetization m z = n+i — 
n_i is conserved, the relevant magnetic energy in an ap- 
plied magnetic field is the second-order (quadratic) Zee- 
man shift of magnitude q = qsB 2 , with B the applied 
magnetic field and qs ~ 277 Hz/G 2 ; The larger (first- 
order) linear Zeeman shift has no influence (it can be 
absorbed in the Lagrange multiplier associated to the 
fixed magnetization). As other spin-changing mecha- 
nisms than collisions are possible, this conservation law is 
only approximate. For example, it no longer holds when 
spin-flips are induced on purpose by applying oscillating 
fields as described above, or for systems with magnetic 
dipole-dipole interactions 



17] . In the absence of such 



applied fields, we find no evidence for violation of this 
conservation law within our experimental limit of a few 
percents. 

We show in Fig.[2]the measured values of no for a range 
of applied magnetic fields B and m z ss 0.4. The popu- 
lation in to = is small at low applied fields and rises 
sharply above a critical value B c before settling at an 
asymptotic value. We have repeated these measurements 
for a wide range of B and to 2 , and generically observed 
this behavior. We show the results in a reconstructed 
contour plot in Fig. [3^- The phase diagram shows un- 
ambiguously the presence of two different phases which 
differ in their spin composition, or more precisely are 
characterized by the absence or presence of condensed 
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FIG. 2: (Color online) Sample data showing the population 
no of the m = Zeeman state versus applied magnetic field 
B, for a magnetization m z ~ 0.4. The solid line is a fit to the 
data using Eq. ((5|. Vertical error bars show statistical uncer- 
tainties on the measured values (one standard deviation). 



atoms in to = 0. 

We now explain the observed behavior of no in terms of 
the competition between the spin- dependent interactions 
and the applied magnetic field (entering quadratically 
through the second order Zeeman effect). The mean- 
field energy functional in the single-mode approximation 
is given by Q 



N 



\S\ -quo- 



(2) 



Here, S = (£\S\Q is the expectation value of the spin 
operator S taken in the normalized spinor £ describing 
the condensate spin wavefunction, and U s denotes the 
spin-spin interaction energy (see Supplementary Mate- 
rial). For antiferromagnetic interactions (U s > 0), no 
applied field (q = 0) and zero magnetization, the spin 1 
BEC realizes a polar, or "spin-nematic" , phase according 
to mean- field theory [l|, Q ■ The spin wave function £ be- 
longs to the family of eigenstates of S ■ n with zero eigen- 
value (and zero average spin), with n a headless vector 
called "nematic director" in analogy with the analogous 
order parameter characterizing nematic liquid crystals. 
When q = 0, any direction n is a possible solution, while 
any positive q favors occupation of the to = state (along 
z) and pins the nematic director in the z direction. 

When m z is non zero, there is a competition between 
the spin-dependent interactions and the quadratic Zee- 
man energy. The constraint of a fixed magnetization is 
essential to understand the spin structure of the conden- 
sate [23| . The BEC spin wavefunction can be parameter- 
ized generically as [ij, [2, 123J 



c = 



no 
no e 



m z ) e^ 1 \ 

'■■4>o 



\ sj ± (1 - n - m z ) e^-i j 



(3) 



We introduced the phases <f> m of the components of C in 
the standard basis. The effect of antiferromagnetic spin- 
dependent interactions (U s > 0) is two-fold: First, they 



lock the relative phase 20o to tt in the minimal 

energy state. Second, they favor the coexistence of the 
to — ±1 component and disfavor mixing them with the 
to = component 0]. As the quadratic Zeeman energy 
favor the latter, the competition between the two results 
in two distinct phases as observed experimentally. 

The equilibrium population no is found by minimiz- 
ing the mean-field energy functional [231 ] . For low q and 
non-zero magnetization m z , spin-dependent interactions 
are dominant, and result in a two-component condensate 
where the Zeeman states m = ±1 are populated (no = 0). 
Following [24| , we will call this phase "antiferromagnetic" 
(AF). When m z — > 0, this gives an "easy-plane" polar 
phase where the nematic director is confined to the x — y 
plane. Above a critical value q c given by 



q c = U a (l- V^^t) 



(4) 



no increases continuously from zero, indicating a second- 
order quantum phase transition. Again following 2J|, we 
call this phase "broken axisymetry" (BA). For large q, the 
energy is minimized by increasing no as much as possible 
given the constraint of a given m z : The spin populations 
therefore tend to n + i = m 2 , no = 1 — m z , n_i = 
for m z > 0. When m z 0, one recovers the easy- axis 
polar phase with all atoms in the m = state along z. 
More generally, the BA state with no 7^ has non-zero 
longitudinal and transverse magnetization (both vanish 
when m z goes to zero), and a nematic director orthogonal 
to the direction of the magnetization vector [25[ . 

We measured the critical line separating the AF and 
BA phases using the following procedure. We bin the 
data according to the measured magnetization, in bins of 
width 0.1 around an average magnetization from m z w 
to m z 1=3 0.8, with residual fluctuations around Sm z s» 
0.02. Each dataset with given magnetization is fitted 
with a function of the form 



no = 



A ,q< q* 



An 



A, q ~ q n > n* 



(5) 



This form ensures the existence of a sharp boundary de- 
termined by q* , a constant background value for low q 
and a well-defined asymptotic value for large q, and re- 
produces the observed data fairly well, as shown in Fig. [5] 
for a specific example with m z w 0.4. At low fields, 
no is not strictly zero but takes values of a few per- 
cents, which can be explained by the presence of a small 
non-condensed fraction (/' « 2 — 3 % per component). 
As such small populations are near our detection limit 
(~ 3 % for the fractional populations, limited by the op- 
tical shot noise associated with the imaging process), we 
do not attempt to determine them and consider in the 
following that the condensate is essentially at zero tem- 
perature. At high fields, no is very close to the expected 
value 1 — m z (see Fig.U^i), again within a few percents. 
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FIG. 3: (Color online) a: Experimental phase diagram showing the population no of the m = Zeeman state versus magne- 
tization m z and applied magnetic field B. The plot shows a contour interpolation through all data points, with magnetization 
ranging from to 0.8. The white line is the predicted critical field B c separating the two phases, deduced from Eq. (|4]) by 
q c = qBBt. b: Theoretical prediction for no at T — K. 



We show in Fig. 0Jd the measured boundary between 
the two phases, which we find in good agreement with the 
prediction of Eq. ^ in the whole range investigated. The 
comparison is made with the value U s /h ~ 65.6 Hz, ob- 
tained from a numerical solution of the Gross-Pitaevskii 
equation using the scattering lengths given in [26| and the 
measured trapping parameters and average atom num- 
ber, and thus does not require any fitting parameter. 
Our results are in line with previous measurements in 
which were restricted to the range m z > 0.5 and 
B > 0.2 G and performed with much larger samples well 
in the Thomas-Fermi regime. Here, we are able to char- 
acterize this transition down to zero magnetization and 
zero applied field, in a system where spin domains (as ob- 
served in [12| during the relaxation towards equilibrium) 
are not expected to form. 

Mean-field theory also quantitatively describes our 
data above the critical line. We compare the calculated 
no directly to the data in Fig. [3^ and b. There is no 
adjustable parameters in this comparison, since the pa- 
rameters used in the theory are either measured or com- 
puted independently. The shape and magnitude of the 
calculated phase diagrams matches well the measured one 
(within 10 % at worst), except very close to the origin 
B sa and m z ss 0. In this corner of the phase diagram, 
we observe larger deviations from the mean field predic- 
tion and correspondingly higher fluctuations in no- We 
will present detailed study on these findings, which go 
beyond the scope of the present paper, in another publi- 
cation. 

In conclusion, we have explored experimentally the 
phase diagram of spin 1 BECs with antiferromagnetic 
interactions. Two phases are found, reflecting the com- 
petition between the spin-dependent interactions and 
the quadratic Zeeman energy. The measurements are 
in quantitative agreement with mean-field theory, which 
quantitatively predicts the phase boundary but also the 
observed spin populations above the transition. One ex- 
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FIG. 4: (Color online) a: Asymptotic value of no for large 
q (determined from Ao + A\ in Eq. [SJ. The solid line shows 
the value 1 — m z expected at zero temperature, b: Mea- 
sured critical field B c versus magnetization. The solid line 
shows the values expected from Eq. (j4|) and q c = qsB\, using 
Us/h ~ 65.6 Hz. The gray area show the uncertainty on the 
theoretical value of B c , dominated by the 15 % uncertainty on 
the spin-dependent scattering length a 3 . For both plots, ver- 
tical error bars show statistical uncertainties on the measured 
values (one standard deviation). 



pects much larger relative changes either at very low tem- 
peratures, where the interplay between spin-dependent 
interactions and quantum depletion are predicted to bc 
important 27, 28|, or conversely at higher temperatures 



where the thermodynamics should be substantially dif- 
ferent from than of the scalar gas Both paths 
provide interesting directions for future work. 
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Supplementary Material 
Sample preparation 

In this section, we give a more comprehensive ac- 
count of the preparation sequence used in the experiment. 
Evaporative cooling is done in two steps as explained in 
[ijjj . starting from a large- volume optical trap that is 
subsequently transferred to a smaller trap with tighter 
confinement (which is used for the experiments described 
in the main text). This sequence allows one to maintain 
a high collision rate throughout the whole evaporation 
ramp. We start from atoms loaded in the large-volume 
trap from a magneto-optical trap (MOT). The loading 
is done at a reduced trap laser power, which was found 
in to optimize the loading. After all MOT lasers 
are switched off, the large- volume trap is compressed by 
ramping up the laser power in 2 s. This increases the col- 
lision rate in the arms of the trap, helps filling the cross- 
ing region and overall provides a better starting point 
for the subsequent evaporative cooling ramp. The laser 
cooling sequence before the compression is found to re- 
sult in a mixed spin state, with spin populations in the 
Zeeman states m = +1,0,-1 in a ratio 0.7 : 0.2 : 0.1, 
approximately. 

To increase the degree of spin polarization, the com- 
pression ramp is done with an additional vertical bias 
field (~ 0.5 G) and magnetic field gradient (20 G/cm). 
As shown in [20(, this results in a spin distillation which 
polarizes the sample into the m = +1 state. The reason 
is that the trapping potential in the vertical direction are 
now slightly different for each Zeeman state, due to the 
potential drop caused by the gradient (see inset of Fig. EJ. 
One can choose a value such that the magnetic potential 
almost compensates gravity for the m = +1 state. The 
772 = state still feels the gravitational potential, and the 
to = —1 state then feels a potential drop twice as strong 
as to = 0. As a result, the effective potential depths for 
to = 0,-1 are slightly reduced compared to the m = +1 
state, and evaporative cooling removes the former atoms 
preferentially. After the spin distillation is complete, we 
obtain a partially polarized cloud with m z ranging from 
fa 0.6 to pa 0.85 depending on the strength of the mag- 
netic field gradient (see Fig. [5k). We found that keeping 
the gradient for longer times was no longer effective to 
increase the polarization further. Our interpretation is 
that as the cloud size becomes too small, the magnetic 
potential drop becomes almost unnoticeable. 

To obtain lesser degrees of polarization, we apply a hor- 
izontal bias field ~ 0.25 G and apply an additional radio- 
frequency (rf) field resonant at the Larmor frequency for 
a variable amount of time. As the atoms move and collide 
in the dipolc trap, their spins quickly decohere, and pro- 
duce a spin-isotropic mixture. By adjusting the strength 
of the rf field, we can adjust the final magnetization at 
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will, as shown in Fig. [5J5. The radio-frequency resonance 
is about 2 kHz wide, presumably limited by inhomoge- 
neous broadening and stray magnetic fields (which are 
estimated to a few mG due to environmental noise). To 
ensure that small frequency drifts do not perturb signif- 
icantly the preparation, the frequency of the oscillating 
field is swept over 20 kHz at a slow rate during the whole 
depolarization sequence. 

After this preparation stage, we perform evaporative 
cooling by reducing the depth of the crossed dipole trap 
until a temperature ~ 10 /iK is reached, at which point 
the cloud is transferred to the final trap with tighter 
focus to boost the spatial density (l9j . This final trap 
is formed by two red-detuned laser beams, one propa- 
gating vertically with a waist (1/e 2 radius) of « 8 /xm 
and the other propagating horizontally with a waist 
w 11 jttm. At the end of the evaporation ramp, where 
the experiments are performed, the trap frequencies are 
Wx, y ,z} = 2tt x (910, 1000,425) Hz. 



Stern-Gerlach expansion 

The populations of the Zeeman states are analyzed af- 
ter expansion in a magnetic field gradient b' = 15 G/cm. 
With an additional bias field B x w 2 G, this produces 
a force along the horizontal x axis that separates the 
m = ±1 clouds from the m = one by a distance 
dsc = ptBVb't 2 /4Afjvai with [1b the Bohr magneton and 
t the expansion time. The factor 7/ takes into account 
the temporal profile of the gradient, which rises in a few 
ms after the beginning of the expansion. Fig. shows 
the vertical trajectory of the atoms, and compares it with 
the one calculated from the measured gradient variations. 
The excellent agreement indicates that the gradient be- 
havior is well understood. After a given expansion time 
(typically £ ~ 3.5 ms), we take an absorption picture of 
the clouds, and count the relative populations. In or- 
der to obtain reliable images, the separation dsc must 
be much larger than the cloud sizes Rt after expansion 
to clearly separate each Zeeman component. In our ex- 
periment, when the trap is switched off instantaneously, 
we typically achieve dsdRt ~ 1 only. This is due to 
the tight trap frequencies, and the resulting fast expan- 
sion. The gradient strength cannot be increased further 
due to technical limitations, and the expansion time is 
also limited by the necessity to keep a sufficiently large 
signal-to-noise ratio to detect atoms in each component. 

We thus resort to a slow opening of the trap, by ramp- 
ing down the laser intensity to approximately l/10th of 
its initial value in 5.5 ms before switching it off abruptly. 
As shown in Fig. |6Jd, this reduces the expansion speed 
(of both the condensate and the thermal gas). At the 
same time, this leaves time for the gradient to settle to 
its maximum value, leading finally to dsa/Rt ~ 10 for 
an expansion time t = 3.5 ms. We have checked that this 



procedure do not affect the measured atom number (see 
Fig. [6t) or condensed fraction. 



Spin interaction energy 

The spin-spin interaction energy (positive for antifcr- 
romagnetic interactions) is given in the SMA by 



U„ = 



Ai:h 2 Na s 



|0(r)| 4 d 3 v, 



(6) 



with niNa the mass of a Sodium atom, a s 0.1 nm the 
spin-dependent scattering length [26| and tfi the single- 
mode wave function. We obtain the latter by solving 
numerically the Gross-Pitaevskii equation [22| 

h 2 - 1 

-^^Acj){r) + -m N ,uj 2 cj) + Ng\cj)\ 2 ( j) = ^. (7) 

We assumed the real, slight anisotropic trap potential 
could be approximated by an isotropic harmonic poten- 
tial, with uj/2tt s» 0.7 kHz the geometric average of the 
three trap frequencies. The spin-independent interaction 
strength is ~g — 4nH a/mssta,, with a 2.79 nm 26[, and in 
accordance with the single-mode assumption we have ne- 
glected spin-dependent interaction terms of order ~ a s /a. 



Conservation of magnetization 

We discuss in this section the key assumption be- 
hind this work, the conservation of longitudinal mag- 
netization. As already discussed, this is true as far as 
short-range interactions are concerned. However, there 
are other weak effects that could in principle relax the 
magnetization. Two main effects come to mind. First, 
a dipole-dipole interaction exists in principle between 
atoms with non-zero spin. These effects are very weak 
compared to short-range spin-dependent interactions. A 
typical relaxation rate due to dipole-dipole interactions 
is less than 0.02 Hz for a fully polarized gas for our pa- 
rameters, assuming the dipolar loss constant is give n by 
the upper bound L 2 ~ 5 ■ 10~ 16 at.cm 2 /s given in [30|. 
Dipolar relaxation can therefore be neglected for the ex- 
periments reported here. Second, the cloud held in the 
optical trap is subject to permanent evaporative cooling, 
leading to a 1/e lifetime around 10 s. If the potential 
depends on the internal state, in general magnetization 
is not conserved (this is the principle behind spin dis- 
tillation [l^]). However, even in this case the remain- 
ing trapped atoms will relax to a new equilibrium state 
through magnetization-conserving collisions, and thus to 
the equilibrium state expected with fixed magnetization. 
In other words, assuming the spin degrees of freedom 
equilibrate faster than the magnetization relaxes, the 
system should adiabatically follow the slow relaxation 
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of magnetization due to evaporation. We note that if 
the potential is spin- independent, and the thermal gas 
isotropic the magnetization will not change on average, 
although one can expect fluctuations to increase in time. 
Experimentally, we found no evidence for these effects, 
which we believe to exist but lie beyond our sensitivity 
(a few percent, limited by the optical shot noise in the 
detection process). 
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FIG. 5: Supplementary Material- (Color online) (a): Spin distillation to prepare samples with high magnetization (m z > 
0.5). The plot shows the magnetization measured for cold clouds, as a function of the magnetic gradient b' in units of 
mNa,g/l-iB, with g the acceleration of gravity and \jlb the Bohr magneton. The inset shows a sketch of the potential energies for 
each Zeeman state along the vertical axis z. The potential drop is exaggerated for clarity, and is smaller than depicted in the 
actual experiment, (b): Depolarization to prepare samples with low magnetizations (m z < 0.5). The time shown corresponds 
to the length of a radio-frequency pulse at the Larmor frequency. 
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FIG. 6: Supplementary Material- (Color online) Free expansion after trap opening, (a): Center-of-mass trajectory along 
the vertical direction z for the m — +1 state; The solid line shows the calculated trajectory taking the measured magnetic 
gradient and gravity into account, (b): Sizes of the expanding clouds after an expansion time for an instantaneous release 
(blue squares) or a smooth release (red circles), (c): Atom number measured for instantaneous (blue squares) and smooth (red 
circles) releases, for various evaporation times. 



